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Abstract 

The dynamics of vortices in a superfluid film at the absolute zero of tempera- 
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pared to a recently proposed quantum-field-theoretic, i.e., "second-quantized" de- 
scription. The theory, consisting of the nonrelativistic effective action of phonons 
linearly coupled to a Chern-Simons term, is used to calculate the one-loop ampli- 
tude for the elastic scattering of phonons from a vortex. 
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1 Introduction 



Vortices play an important role in several two-dimensional condensed matter systems 
at low temperatures such as superfluids and superconductors, fractional quantized Hall 
systems [[[J and Josephson junction arrays [|J|. They are also believed to play an im- 
portant role in the quantum phase transitions these systems undergo |], Bj. Two- 
dimensional vortices can ideally be represented as pointlike and are, close to the ab- 
solute zero of temperature, governed by the laws of quantum physics. The quantum 
mechanics of such pointlike vortices, which is known since long (see, for example, 
Ref. [^]), is unusual. This is because classically, instead of obeying Newton's equa- 
tions of motion, the coordinates of a vortex center obey Kirchoff 's equations of motion 
[ph. As will be recalled in Sec. ^, this is connected to a peculiar canonical structure 
known to arise also in the problem of a charged particle confined to move in a plane 
subject to an external magnetic field [ph. As a result, the classical dynamics of a vortex 
in a two-dimensional fluid resembles that of a charged particle in a plane under the 
influence of a magnetic field [Q] . 

In recent years, various quantum-mechanical formulations of vortex dynamics in 
superfluid films have been put forward which use this analogy [|[ [To| . In these 
so-called dual approaches (see Sec. [|), vortices are still described by delta functions 
representing their worldlines, but their interaction is now pictured as being mediated 
by a photon — the quantum of the electromagnetic field — and not, as in the original 
formulation, by a phonon — the quantum of the Goldstone field associated with the 
spontaneously broken U(l) symmetry of a superfluid. The analogy can be understood 
| [l2| | by comparing the continuity equation 

d n + V-j = 0, (1) 

where n is the particle number density and j the particle number current of the super- 
fluid, to the Maxwell equation (Gaussian units will be used throughout): 

-d H + V x E = 0. (2) 

c 

It is to be noted that in 2+1 dimensions, the magnetic field H is a scalar, and VxE = 
eijdiEj, with ey = 1,2) the antisymmetric Levi-Civita symbol (ei2 — 1). It 
follows that with the replacements 

n^iJ/$ , ji^-reijEj, (3) 

where <l>o = hc/e is the magnetic flux quantum, the continuity equation goes over 
into the Maxwell equation. The fundamental constants are included in ([|) to match the 
dimensions on both sides of the replacements. Since the integral $(5) = J s d 2 xH 
represents the magnetic flux penetrating the surface S, the ratio Hj <£>o is the flux num- 
ber density n®. In other words, the duality transformation maps the particles of the 
original formulation onto flux quanta in the dual formulation. 

In this paper, we take an entirely different approach. Instead of considering a 
quantum-mechanical, i.e., "first-quantized" description of vortices in a superfluid film 
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at the absolute zero of temperature, we study a quantum-field-theoretic, i.e., "second- 
quantized" description recently put forward by these authors [|l3]] . Vortices are here not 
described by delta functions representing their worldlines as in the quantum-mechanical 
approach, but by a nonsingular quantum field. The quantum field theory fulfills two 
stringent constraints. The first one arises because at zero temperature the entire system 
is superfluid, and can, therefore, be considered an ideal fluid. According to a theorem 
due to Helmholtz [pp, a vortex in an ideal fluid moves with the fluid, and, hence, has no 
independent dynamics. The second constraint stems from the following observation. 
An electron traversing a closed path in a magnetic field accumulates an Aharonov- 
Bohm phase given by (2ir times) the number of flux quanta encircled [Q. The above 
analogy then asserts that an elementary vortex traversing a closed path in a superfluid 
film also accumulates a geometrical phase given by (2ir times) the number of particles 
encircled [|l5|]. It has been shown in [|l3|| (see also Sec. that both these constraints 
are satisfied by a specific Chern-Simons theory. Two of the main characteristics of a 
Chern-Simons theory in general are that, first, it contains a vector field which has no 
independent dynamics, and that, second, it imparts flux to particles. By representing 
vortices by a Chern-Simons field, we assured that they have no independent dynamics, 
and that they see the particles as sources of geometrical phase. To connect with the 
standard quantum-mechanical approach, we introduce in Sec. ^external vortices in the 
field theory and show that, in this way, the known results are reproduced. 

The specific Chern-Simons theory proposed in [jjjj] is based on the effective theory 
of phonons describing a superfluid film at low energies and small momenta (see Sec. 

. The effective Lagrangian is the most general one consistent with the symmetries 
of the problem, in particular Galilei invariance and invariance under global U(l) trans- 
formations. In Sec. U we use this effective theory extended by the Chern-Simons part, 
representing the vortices, to investigate one-loop contributions to the elastic scattering 
of phonons by a vortex. We apply dimensional analysis to identify the relevant Feyn- 
man graphs. To the order in which we are working, we do not encounter ultraviolet 
divergences in the effective theory; only a logarithmic infrared divergence arises. This 
introduces an arbitrary scale factor into the problem. Often such a factor appears in- 
stead when eliminating ultraviolet divergences by a redefinition of certain parameters 
of the effective theory. The one-loop contributions evaluated in this paper give the first 
correction to an old result for the scattering of a sound wave from a vortex obtained by 
Pitaevskii in the Born approximation [Q. 

Notation In the main text we will keep track of the fundamental constants. We 
will follow standard notation and use p as momentum variable and k as wave-vector 
variable; they are related via p = Kk, where h = h/2ir with h Planck's constant. 
The (angular) frequency will be denoted by u and the corresponding energy by p^; 
Po = Tilo. A spacetime point will be indicated as a; = — (t, x), /i = 0,1, ... ,d, 
with d the number of space dimensions, the energy-momentum as p = p„ = (p , p), 
and k = k^ = (uj, k). 
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2 Effective Theory for Sound Waves 



In this section, we briefly recall the nonrelativistic effective action describing sound 
waves at low energy and small momentum in two space dimensions [|l7]]. As starting 
point to describe the system, we take the microscopic model | 



C = <j>* [ihd - e(-iV) + - A o |0| 4 , (4) 

where the complex scalar field 4>{x) describes the atoms of mass m constituting the 
liquid, e(—ihV) = — ?i 2 V 2 /2m is the kinetic energy operator, and /io is the chemi- 
cal potential. The last term with a positive coupling constant, Ao > 0, represents a 
repulsive contact interaction. 

The Lagrangian has a global U(l) symmetry, under which the matter field trans- 
forms as 

4>{x) -> <t>'{x) = e ia 4>(x), (5) 

with a the transformation parameter. At zero temperature, this global symmetry is 
spontaneously broken by a nontrivial ground state. This can be easily seen by consid- 
ering the potential energy 

V = — a^oI^I 2 + A |^| 4 , (6) 
which has, for /iq > 0, a minimum away from the origin <fi = at 

\4>\ 2 = (7) 

Z Ao 

Since the total particle number density n(x) is given by 

n(x) = \cp{x)\\ (8) 

the quantity uq := \<J>\ 2 physically represents the number density of particles contained 
in the condensate. 

To account for the nontrivial ground state, we shift <f> by the (complex) constant <f> 
and write ^ 

cj>(x)=e iv ^[4> + j>{x)}. (9) 

The scalar field (f(x) is a background field representing the Goldstone mode of the 
spontaneously broken global U( 1 ) symmetry. In terms of the new variables, the quadratic 
terms of the Lagrangian (Q) may be cast in the matrix form 

Co = ±&M $, l=fU, (10) 

with 



, ihdo - e + ^0 - U - 4A O |0| 2 -2A ^ , 

Ma ~ { -2A O 0* 2 -ihd -e + fio-U-4X \ Ia 1 ' (U) 



and U is the Galilei-invariant combination [| 

1 



U = hd„<p+ ~-{KS7<p) 2 . (12) 
2m 
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In writing (|TT]), we have omitted a term of the form V 2 ip containing two derivatives 
which is irrelevant in the regime of low momentum in which we are interested. We 
have also omitted a term of the form V</? ■ j, where j is the Noether current associated 
with the global U(l) symmetry, 

j=2^(tf*W-W«. d3) 

This term, which after a partial integration becomes —ipV ■ j, is also irrelevant at low 
energy and small momentum because in a first approximation the particle number den- 
sity is constant, so that the classical current satisfies the condition 

V-j = 0. (14) 

The spectrum E(p) obtained from the matrix Mq with the background field U set to 
zero is the famous single-particle Bogoliubov spectrum jl9|], 

E(p) = v / e 2( p ) + 2 / i e(p) 

= ^e 2 (p)+4A o |0|2 e (p). (15) 

The most notable feature of this spectrum is that it is gapless, behaving for small mo- 
mentum as 

£(p)~uo|p|, (16) 

with uq = y/na/m a velocity which is sometimes referred to as the microscopic sound 
velocity. It was first shown by Beliaev Jl8| ] that (in three dimensions) the gaplessness 
of the single-particle spectrum persists at the one-loop order. This was subsequently 
proven to hold to all orders in perturbation theory by Hugenholtz and Pines [pp[]. For 
large momentum, the Bogoliubov spectrum takes a form 

£(p)~e(p) + 2A o |0| 2 (17) 

typical for a nonrelativistic particle with mass m moving in a medium. To highlight 
the condensate, we have chosen here the second equality in (|l^), where /j,q is replaced 
with2A H 2 . 

The gaplessness of the single-particle spectrum ( fl"5] ) is a result of Goldstone's theo- 
rem. This can, for example, be seen by considering the relativistic version of the theory. 
There, one finds two spectra: one corresponding to a massive Higgs particle which in 
the nonrelativistic limit becomes extremely heavy and decouples from the theory, and 
one corresponding to the Goldstone mode of the spontaneously broken global U(l) 
symmetry [pT||. The latter reduces in the nonrelativistic limit to the Bogoliubov spec- 
trum. Another way to see this is to couple the theory (Q) to an electromagnetic field. 
The single-particle spectrum then acquires an energy gap. This is what one expects to 
happen to the spectrum of a Goldstone mode when the Higgs mechanism is operating. 
In other words, the single-particle spectrum ( p"5| ) is identical to that of the Goldstone 
mode which physically represents the collective density fluctuations in the superfluid. 
The equivalence of the single-particle excitation and the collective density fluctuation 
has been proven to all orders in perturbation by Gavoret and Nozieres p2[|. 
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Figure 1: Graphical representation of the effective theory f^j. The symbols are ex- 
plained in the text. 



To find the effective theory governing the superfluid at low energy and small mo- 
mentum from (|l0|), we integrate out the fluctuating field $ Jl7[], The first two terms of 
the effective theory are graphically represented by Fig. [ij A line with a shaded bubble 
inserted stands for ifi times the full Green function G and the black bubble denotes i/h 
times the full interaction T of the $-field with the background field U which is denoted 
by a wiggly line. Both G and T are 2x2 matrices. The full interaction is obtained from 
the inverse Green function by differentiation with respect to the chemical potential, 

r--?f-. as) 

dpi 

This follows because U, as defined in (|l2|), appears in the theory only in the combi- 
nation [I — U. To lowest order, the inverse Green function is given by the matrix Mq 



in ( 1 1 ), so that the vertex describing the interaction between the $ and [/-fields is i/h 
times minus the unit matrix. Because in terms of the full Green function G, the particle 
number density reads 

ih fd d+1 p 



6 tr J jzrt G (p)> ( 19 ) 

where d d+1 = dpod d p, it follows that the first diagram in Fig. |l| stands for —nil. 
The bar over n is to show that the particle number density obtained in this way is a 
constant, representing the density of the uniform system with U(x] set to zero. The 
second diagram without the wiggly lines denotes i/h times the (0 0)-component of the 
full polarization tensor, IIoo, at zero energy transfer and low momentum q, 

|limnoo(0,q)=-ilimtr J ^—^-GTG (po,p + q). (20) 
The factor i is a symmetry factor which arises because the two Bose lines are identical. 



We proceed by invoking an argument due to Gavoret and Nozieres [ 22 1 to relate the 
left-hand side of ( j20| ) to the sound velocity. By virtue of relation ( |l8| ) between the full 
Green function G and the full interaction T, the (0 0)-component of the polarization 
tensor can be cast in the form 



ih f d o 9G 

limIT 00 (0,q) = -Yl^try^rG — G(p ,P + q) 

ih d r [d d+1 p„, 

= y^^o tr J /^ G(po ' p+q) 
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On 
Oft, 



i an 



(21) 



where f2 is the thermodynamic potential and V the volume of the system. The right- 
hand side of ( p"i| ) is n 2 n, with k the compressibility. Because it is related to the macro- 
scopic sound velocity c via 

«=— ^, (22) 
mnc z 

we conclude that the (0 0)-component of the full polarization tensor satisfies the so- 
called compressibility sum rule of statistical physics 



lim n oo (0,q) 



(23) 



Putting the pieces together, we infer that the diagrams in Fig. [j] stand for the effective 
theory 



C e fi 



hd <p+—(h\7ip) 2 
Im 



2mc 2 



hd Q ip + —(h\7<p) 2 
Am 



(24) 



where we recall that n is the particle number density of the fluid at rest. This nonlinear 
theory describes a nonrelativistic sound wave, with the dimensionless real scalar field cp 
representing the Goldstone mode of the spontaneously broken global U(l) symmetry. It 
has the gapless dispersion relation w 2 (k) = c 2 k 2 , where to is the (angular) frequency, 
k the wave vector, and c the sound velocity. This effective theory gives a complete 
description of the superfluid valid at low energies and small momenta. There are of 
course higher-order terms, but they need to be included only as one goes to higher 
energies and momenta. The same effective theory we discussed here appears also in 
the context of (neutral) superconductors and that of classical hydrodynamics [f23|]. 

The form of the effective theory (24) can also be derived from general symmetry 
arguments [24 25j|. The basic idea is that the presence of a gapless Goldstone has to 
be reconciled with Galilei invariance, which demands that the mass current and the 
momentum density are equal. This leads to the conclusion that the U(l) Goldstone 
field <p can only appear in the combination (|l2|). To obtain the required linear spectrum 
for the Goldstone mode it is necessary then to have the form (p4[). 

The particle number current that follows from (24) reads 



n(x) 



^-\%dw{x) + ^-[KV<p{x)] 



mc 
n(x)v(x), 



2m 



(25) 
(26) 



where v = (h/m)Vf is the superfluid velocity field. Equation (J25[) reflects Bernoulli's 
principle which states that in regions of rapid flow, the density and therefore the pres- 
sure is low. 

The diagrams of Fig. [l] can also be evaluated in a loop expansion to obtain explicit 
expressions for the particle number density n and the sound velocity c to a given or- 
der jl7|]. In doing so, one encounters — apart from ultraviolet divergences which can 
be eliminated by renormalization — also infrared divergences because the Bogoliubov 
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spectrum is gapless. When, however, all one-loop contributions are added together, 
these divergences are seen to cancel [p"7|]. One finds for d = 2 to the one-loop order 

n = % c*=»=2±n, (27) 
2 A mm 

where A and fj, are the renormalized coupling constants (see below). The second ex- 
pression in (|7[) for the sound velocity is appropriate when, as is often the case in ex- 
periment, the particle number is fixed. The first one, featuring the chemical potential, 
is appropriate in the presence of a reservoir with which the system can freely exchange 
particles, so that only the average particle number is fixed. 

Given the form of the effective theory, the particle number density and sound ve- 
locity can be more easily obtained directly from the thermodynamic potential 57 via 

i an i i md 2 n 

V o/i <r V n d[L A 

where V is the volume of the system. In this approach, one only has to calculate 
the thermodynamic potential which at zero temperature and in the approximation we 
are working is given by the effective potential V c s corresponding to the theory (^): 
fi = J d d a;Vcff. To the one-loop order, the effective potential for a uniform system 
reads 



with E(p) the gapless Bogoliubov spectrum (jl5|). The integral over the loop momen- 
tum in arbitrary space dimension d yields 

v Mo r(l-rf/2)r(rf/2 + l/2) d/2 d/2+1 

V ° ff = ~TXo - 2.(^/^/2 + 2)™ M ° ' (30) 

where we used the integral representation of the Gamma function 

oo 

1 Wdr (31) 



a z T(z) J t 
o 

together with dimensional regularization to suppress irrelevant ultraviolet divergences. 
By this we mean contributions which, for regularization with a momentum cutoff A, 
would diverge with a strictly positive power of A. Note that we, following Ref. [p6|, 
first carried out the integrals over the loop energies, and then analytically continued the 
remaining integrals over the loop momenta to arbitrary space dimensions d. 
Expanded around d = 2, Eq. ( |30| ) gives 

V eS = -^--^^ + 0(e% (32) 
4A AtzTi e « 

where e = 2 — d is the deviation from the upper critical dimension d = 2, and n is 
an arbitrary renormalization group scale parameter, with the dimension of an inverse 
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length. The right-hand side of ( p2| ) is seen to diverge in the limit d — ► 2. The theory 
can be rendered ultraviolet finite by introducing a renormalized coupling constant A via 



1 1/1 



-2 

■ e 



(33) 



where A = A/k £ . Its definition is such that for arbitrary d, A has the same engineering 
dimension as Ao in the upper critical dimension d — 2. As renormalization prescription 
we used the modified minimal subtraction. The beta function /3(A) follows as [E7fl 



/3(A) = k ^ 

OK 



= -eA+-^A 2 . (34) 



In the upper critical dimension, this yields only one fixed point, viz. the infrared-stable 
fixed point A* = 0. Below d = 2, this point is shifted to A* = enh 2 jm. 

It follows from ( f32| ) that the chemical potential is not renormalized to this order. 
Incidentally, from the view point of renormalization, the mass m is an irrelevant pa- 
rameter in nonrelativistic theories and can be scaled away. 

The most remarkable aspect of the effective theory (Q) is that it is nonlinear. The 
nonlinearity is necessary to provide a Galilei-invariant description of a gapless mode in 
a nonrelativistic system. Since the Goldstone field in (Q) is always accompanied by a 
derivative, we see that the nonlinear terms carry additional factors of h\k\/mc, with |k| 
the wave number. They can therefore be ignored provided the wave number is smaller 
than the inverse coherence length £ = %/ mc, 

|k| < 1/e (35) 

For example, for 4 He the coherence length, or Compton wavelength, is about 10 nm. 
In this system, the bound ([[5]), below which the nonlinear terms can be neglected, 
coincide with the region where the spectrum is linear and the description in terms of a 
sound mode is applicable. 

We close this section discussing the apparent mismatch in the number of degrees of 
freedom in the normal and superfluid state. Whereas the normal state is described by a 
complex field <fr, the superfluid state is described by just a single real scalar field ip. The 
resolution of this paradox lies in the spectrum of the modes [f28|]. In the normal state, 
the spectrum E(p) = p 2 /2m is linear in E, so that only positive energies appear in the 
Fourier decomposition of the field cf>. One needs therefore a complex field to describe 
a single particle, as is well known from standard quantum mechanics. In the superfluid 
state, where the spectrum E 2 (p) = c 2 p 2 , is quadratic in E, the Fourier decomposition 
of the field <p contains positive as well as negative energies. As a result, a single real 
field suffices to describe this mode. In other words, although the number of fields is 
different, the number of degrees of freedom is the same in the normal and superfluid 
state. 



3 Quantum Mechanics of Vortices 

In this section, we discuss the conventional quantum-mechanical, or "first-quantized" 
description of vortices. We will be working at the absolute zero of temperature, so that 
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the entire liquid is superfluid. In the microscopic theory (Q), the asymptotic solution of 
a static vortex with winding number w centered at the origin is well known | 

A* ft t-2 w2 \ Awe , //-i f 1 



fx^ 1 -***)* +0 {*)> (36) 

where 6 is the azimuthal angle and £ is the coherence length introduced above ( |35| ) 
which because of ( p7[ ) can also be written as £ = h/y/rnjl. The density profile n(x) in 
the presence of this vortex follows from taking |0(x)| 2 . 

Let us now discuss how vortices can be incorporated in the effective theory. To this 
end we follow Kleinert [ ^o| ] and introduce a so-called plastic field ip F — (tp P , <£> p ) in 
the effective theory (24) via minimally coupling to the Goldstone field: 

d^^d^ + ip*, (37) 

with <9 M = [do, —V). If there are N vortices with winding number w a (a — 1. • • • , N) 
and centered at X 1 (i), ■ ■ • , X JV (t), the plastic field satisfies the relation 

V x (p F (x) = -2ir^2w a S[x~X. a {t)}, (38) 

a 

so that we obtain for the velocity field 

Vxv = ^ 7a %-X a (f)], (39) 

OL 

as required. Here, j a — (h/m)w a , with w a — 0, ±1, ±2, . . ., denotes the circulation 
of the ath vortex which is quantized in units of h/m. A summation over vortex labels 
will always be made explicit. The combination d^(p + (p F is invariant under the local 
gauge transformation 

tp(x) -> tp(x) + a{x); <p F (p F - df,a(x), (40) 

with ip F playing the role of a gauge field. 

In the gauge </?q = 0, a solution of Eq. ( ^8| ) is given by 

<Pi( x ) = Zireij^WcSjix, L a (t)\, (41) 

a 

where is the antisymmetric Levi-Civita symbol in two dimensions, with €12 = 1, 
and S[x, L a (t)] is a delta function on the line L a (t) starting at the position X Q (£) of 
the ath vortex and running to spatial infinity along an arbitrary path: 

Si[x,L a (t)] = f d W <y(x-y). (42) 

JL a (t) 

Let us for the moment concentrate on static vortices. The field equation obtained 
from the effective theory ( [24] ) with Vip replaced by the covariant derivative Vip — ip F 
and doip set to zero simply reads 

V • v = 0, or V ■ (Vtp - <p p ) = 0, (43) 
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when the fourth-order term is neglected. It can be easily solved to yield 

V(x) = - /d 2 yG(x-y)V V(y), (44) 
where G(x) is the Green function of the Laplace operator in two space dimensions 

/(jSi ikx -i 

i.e., V 2 G(x) = — £(x). For the velocity field we obtain in this way the well-known 
expression [Q] 

«i(x) = e l3 Y, 7«0j f d 2 yG(x - y)% - X Q (i)] 

a 1 1 

which is valid for x sufficiently far away from the vortex cores. Let us now specialize 
to the case of a single static vortex centered at the origin. On substituting the corre- 
sponding solution in (|25|), we find for the density profile in the presence of a static 
vortex asymptotically 

n(x)=nfl-£ 2 ^V (47) 



2x 2 , 

Since this density profile is identical to the one obtained in the microscopic theory [see 



Eq. (pq)], this exemplifies that, with the aid of plastic fields pOfl , vortices are correctly 
accounted for in the effective theory. 

Let us proceed to investigate the dynamics of vortices in this formalism and derive 
the action governing it. We consider only the first part of the effective theory (|4|). In 
ignoring the higher-order terms, we approximate the superfluid by an incompressible 
fluid for which the particle number density is constant, n{x) — n, see Eq. (^5|). This is 
an approximation often used. We again work in the gauge ip^ = and replace Vy> by 
the covariant derivative Vip — y p , with the plastic field given by (f38|). The solution of 
the resulting field equation for ip is again of the form (fw)), but now it is time-dependent 
because the plastic field is. Substituting this in the action S e s — J dtd 2 xC c s, we find 
after some straightforward calculus 



S e g = mn / dt 



2tt 

Q</3 



(48) 



where a dot over a symbol denotes the time derivative. This action yields the well- 
known equations of motion for point vortices in an incompressible two-dimensional 
superfluid [00]: 

x * {t) = ~2^ €ij £ la |x/3(;)-x«(i)| 2 • (49) 
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Note that X?(t) = v t [X^(t)], where v(x) is the superfluid velocity (^) with the 
time-dependence of the vortex centers included. This nicely illustrates the result due 
to Helmholtz for ideal fluids [|7|] we alluded to in the Introduction, that a vortex moves 
with the fluid, i.e., with the local velocity produced by the other vortices in the system. 
Experimental support for this conclusion has been reported in Ref . [[32]] . 

The action (^8j) leads to a twisted canonical structure. To display it, let us rewrite 
the first term of the corresponding Lagrangian as 

L x = -mnJ2^X^X^, (50) 

a 

where we ignored a total derivative. It follows that the canonical conjugate to the 
second component of the center coordinate X Q is its first component [|| 

dL 

-± = -mn 7 „Xf. (51) 

OA 2 

This implies that phase space coincides with real space and it gives the commutation 
relation 

\X?,X5\ = — £ 2 S a(} , (52) 

W a 

where 

^=1/V2^ (53) 

is a characteristic length. Its definition is such that 2tt£ 2 is the average area occupied 
by a particle of the superfluid film. The commutation relation leads to an uncertainty 
in the location of the vortex centers 

AX i AX £ > (54) 

which is inverse proportional to the particle number density. 

Elementary quantum mechanics | |3~3| ] tells us that in the quasi-classical approxima- 
tion to each unit cell (of area h) in phase space there corresponds one quantum state. 
That is, the total number of states is given by 

# states = i J dp dq, (55) 

where p and q are a pair of canonically conjugate variables, and the integral is over the 
entire phase space. For a two-dimensional quantum vortex with winding number wp, 
this implies that the number of states it can be in is 

# states = \wp \ nS, (56) 

where S is the surface area of the sample. In other words, every particle in the super- 
fluid film makes (wp times) an additional state available to the vortex. 

This phenomenon that phase space coincides with real space is known to also arise 
in the Landau problem of a charged particle confined to move in a plane subject to a 
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Figure 2: Duality transformation. 



(constant) background magnetic field. There, it leads to the well-known result that the 
number of states available to the charged particle in each Landau level is 

TJ 

# states = \eg\—S, (57) 
he 

where H is the magnetic field component perpendicular to the plane, and 

ep = vpe , (58) 

(vp = 0, ±1, ±2, . . .) is the electric charge of the particle given as a multiple of the unit 
of charge eo(> 0). In terms of the magnetic flux quantum $o = hc/eo, the number of 
states can be rewritten as 

TJ 

# states = \vp\—S = \vp\n®S, (59) 

where fi® is the flux number density. Hence, whereas the number of states for vortices 
in a superfluid film is determined by the particle number, here it is determined by the 
flux number. This agrees with the replacement (|j) discussed in the Introduction. 

As remarked there, the analogy between the two problems is the basis of a much 
used duality transformation. If we collectively refer to the electric charges in the Lan- 
dau problem and the particles in the superfluid as 'charges', and to the vortices in the 
superfluid and the flux quanta in the Landau problem as 'vortices', then the duality 
transformation interchanges charges and vortices (see Fig. fy. Using this analogy, we 



see that the characteristic length (53) translates into 



l B = 1/V27TO®, (60) 

which is precisely the magnetic length of the Landau problem. 

We next turn to the geometrical phase In the two-dimensional Landau prob- 
lem, when a charged particle, say the /3th, is moved adiabatically around a close path 
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Tp, its wavefunction accumulates an extra Aharonov-Bohm phase factor 7(1^3) given 
by the Wilson loop: 



W(Tp) = exp [t7(I»] = exp 




exp 



2nivp 



HS(Tp) 
$0 



(61) 



where A is the vector potential describing the background magnetic field and HS(Tp) 
the magnetic flux through the area 5 f (r^) spanned by the loop Tp. The geometrical 
phase 7(r i a) in ( [sT] ) is seen to be (2-nvp times) the number of flux quanta enclosed by 
the path Tp. 

Because of the above analogy, it follows that the geometrical phase acquired by the 
wavefunction of a vortex when it is moved adiabatically around a closed path in the 
superfiuid film is (2Ttwp times) the number of particles enclosed by the path [O]. 

We have up to now considered a constant magnetic background field and an in- 
compressible superfiuid characterized by a constant particle number density. A more 
general analysis of the geometrical phase for a compressible superfiuid was carried 
out by Haldane and Wu [15]. Their starting point was an Ansatz for the multivortex 



wavefunction of the interacting Bose condensate 



(62) 



where ^° describes the condensate in the absence of vortices. The operator A(X) 
introduces an elementary vortex centered at X, 



A(X)=l[[(x a 1 -X 1 )+i(x a 2 -X 2 )], 



(63) 



where x a denotes the coordinate of the ath particle, and the factor /(X) accounts 
for the change in the particle number density owing to the presence of this vortex. 
Repeating steps first carried out by Arovas, Schrieffer, and Wilczek [[36] in the context 
of the fractional quantized Hall effect, they cast the geometrical phase [33], 



1 (Tp) = ii dX^ 3 - <tf|V X ftf) 
Jr ff 



(64) 



picked up by the wavefunction when the /3th vortex is moved along a closed path Tp 
in the form 

7(I» = Jd 2 x<£ dX p • V x [inflx-X"!)] n(x;X a ,X /3 ). (65) 

Here, n(x; X", X' 3 ) is the particle number density in the presence of the transported 
vortex as well as the other vortices centered at X" (a (3). Writing 



i(x; X Q , X ) = n(x; X Q ) + Sn(x - X' 3 ), 



(66) 
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where n(x;X Q ) is the particle number density in the absence of the vortex being 
moved, they concluded that, apart from corrections owing to residual vortex interac- 
tions that become small in the dilute-vortex limit, the geometrical phase is given by 



7(I» =2tt / d 2 xn(x;X a ). 



(67) 



That is, for a compressible superfiuid a similar result is found as in the simpler case of 
an incompressible superfiuid, in that the geometrical phase acquired by an elementary 
vortex being moved around a closed path is (2tt times) the average number of particles 
enclosed by the path [|T^]. 

The last topic we wish to discuss in this section on the quantum mechanics of 
vortices is the elastic scattering of a sound wave from a vortex. The subject, first 
studied by Pitaevskii Jl6[|, has recently received considerable attention (see Ref. p7| ] 
and references therein). To obtain the hydrodynamic equation governing this scattering 
process, we first record the nonlinear field equation obtained from the effective theory 
(H© 

<9o/i + c 2 V • v = <9 v 2 + iv ■ Vv 2 , (68) 

with /i = —jh/ m)do(p the chemical potential per unit mass. If we ignore the nonlinear 
terms, Eq. ( pq ) becomes the more familiar wave equation 

dl<p - c 2 V 2 tp = 0. (69) 

Let us next introduce a freely moving, i.e., non-pinned vortex. Driven by the sound 
wave, the vortex will oscillate around some point, x = say. The velocity field in the 
presence of a moving vortex is obtained from the static solution [see Eq. ([46|)1 

"iW = (70) 

Z7T X 

by replacing the coordinate x with x — v^i, where Vi(i) is the velocity of the vortex. 
This implies that 

<9 v v (x- v L t) = -v L ■ Vv v (x- v L t), (71) 

where v is given an index v to show that it is the velocity field produced by the vortex. 
Since the solution v v is curl-free outside the vortex core, the right-hand side of Eq. 
( |7l| ) may be written as — V(v£ ■ v v ) there. The hydrodynamic equation governing the 
elastic scattering of a sound wave from the free vortex now follows from writing the 
velocity field as 

v(x) = v v (x- w L t) + — V<f(x), (72) 
m 

with (p describing small variations around the oscillating vortex solution, and noting 
that (^TJ) requires that we write for the chemical potential per unit mass 

K x ) = v £ ' v v(x- Vx,t) - — do<p(x). (73) 

m 



Substituting this in the field equation (pq), we obtain to linear order in (p the hydrody- 
namic equation [^7|[, 

d&(x) - c 2 V 2 0(x) = -v v (x) • Vd [2ip(x) - (p(t, 0)], (74) 
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where we approximated v v (x — vj,i) by v v (x). To linear order in (p this is allowed 
since the vortex, being driven by the sound wave, has a velocity 

v L (t) = — V£(t, x = w L t) « —V<p(t, 0). (75) 
m to 

We also neglected a term quadratic in v v which is justified because |v v | << c out- 
side the vortex core [ 16 1 . The first term at the right-hand side of ( f74j ) stems from the 



nonlinear term <9ov 2 in the field equation ( |58[) 

Equation (|74|) has been used by Sonin [ [37| ] as a basis to study the elastic scattering 
of a sound wave from a free moving vortex. In the Born approximation, he found 
as two-dimensional scattering amplitude f(9, |k|) the result first derived by Pitaevskii 

m|k|) = iJMAe^^|, (76) 

2 V Z7T TOC 1 — cos 

where 9 is the angle between the incoming and the scattered sound wave, and |k| their 
wave number. 

In the next section, we discuss a quantum field theory of vortices recently proposed 
by these authors ijjjj] and show that it reproduces the results discussed here. 



4 Quantum Field Theory of Vortices 

In the previous section we have seen that vortices in a two-dimensional superfluid at 
the absolute zero of temperature can be considered as point particles which are subject 
to the laws of quantum physics. It is therefore natural to ask whether vortices also 
admit a quantum-field-theoretic, or "second-quantized" description. In a recent letter, 
we have argued they indeed do and proposed the following quantum field theory [OR: 



C = £ c ff - ena + -j ■ a + £ cs . (77) 

c 

It consists of the nonrelativistic effective Lagrangian ( p4| ) describing the superfluid 
without vortices, linearly coupled via the particle number current = (n, j) to a 
vector field = (ao, a) describing the vortices. This field is governed solely by a 
Chern-Simons term 

£ C s = 7r a x d o a - fl o v x a > (78) 
2c 

with the charge e determined by [ pj[ ] 

e 2 = he. (79) 

Our convention is such that, in contrast to j^, the dimensions of the time and space 
components of the vector field a M are the same. It is important to note that the coupling 
of the sound mode to the Chern-Simons field is linear and not minimal as in gauge the- 
ories. The coupled theory ([77]) therefore does not possess a gauge invariance involving 
a simultaneous local gauge transformation of the Chern-Simons field and the matter 
field. To appreciate why a minimal coupling is not feasible here, note that this would 



15 



inevitably result in the disappearance of the gapless sound mode because of the Higgs 
mechanism. This would disagree with the physics we wish to describe, namely the 
dynamics of vortices in a compressible superfiuid film. For the gapless sound mode to 
survive a coupling to a gauge field, the coupling cannot be minimal so that there is no 
local gauge invariance in the full theory. 

The appearance of a Chern-Simons term in ([77]) is expected for two reasons. First, 
as we discussed in the previous section, vortices in a superfiuid at the absolute zero 
of temperature always move with the fluid. In other words, they have no independent 
dynamics. This is precisely also a basic property of a Chern-Simons term. It comes 
about because the zeroth component ao of the vector field is a Lagrange multiplier, 
demanding the equality 

V x a = -en, (80) 

or when integrated 

$ = -eN, (81) 

where N = J d 2 x n is the particle number and $ = J d 2 xV x a the "magnetic" flux 
associated with the Chern-Simons field. In the Coulomb gauge V • a = 0, Eq. ( ^o| ) can 
be easily solved to yield 



a(a;) = -eVx J d 2 yG(x - y)n(t, y), 



(82) 



where the Green function G(x) is given in (p3[). The solution shows that the Chern- 
Simons field is entirely determined by the particle number density n and that it has, 
therefore, no independent dynamics. 

The second reason why a Chern-Simons term is expected to appear is that it en- 
codes the geometrical phase acquired by a vortex when it winds around a boson. The 
quantum-mechanical analog of such a term, representing the linking number [|8]] of 
a closed boson and vortex trajectory, was introduced in the problem by Arovas and 
Freire JTT{]. To calculate the geometrical phase j(T) we again evaluate the Wilson loop 
VK(r) = exp[i7(r)] obtained by integrating the Chern-Simons field a around a closed 
path r in the superfiuid, 

7 (D = ± I dl • a (as) = & — [ d 2 y <j> dl • [V x x lnflx - y\)n(t, y)], (83) 
he J r he J J r 



where we substituted the explicit form ( |82| ) for a. This geometrical phase is pre- 
cisely the one in obtained by Haldane and Wu Jl^]. Whereas in their quantum- 
mechanical description, the vortex being taken around the closed path V sees the encir- 
cled bosons as sources of geometric phase, here this counting is provided by the flux 
imparted to a particle by the Chern-Simons term [see Eq. d8lj)]. 

We next introduce external point vortices into the theory to determine their action 
and compare it with the action ( |48| ) obtained from the quantum-mechanical description 
of vortices. As we did there, we approximate the superfiuid by an incompressible fluid 
and consider only the following terms of the Lagrangian ([77|): 



nd tp + -^—{nvifi) 2 

2m 



-j • a oxt . (84) 

c 
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The particle number current j reads j = nv in this approximation, and the Chern- 
Simons field a cxt is given by (K3) with n replaced by the external vortex density 



"cxt 



(x) = ^> Q <5[x-X Q (i)]. 



(85) 



Explicitly, 



or 



— > w a v arctan „^ . . 

2n ^ \xi-Xf(t) 

-*?(*) 



(86) 



(87) 



where X" (t) denotes the center of the ath vortex with winding number w a . Note that 
the gauge V • a = is satisfied by the solution ( |87| ) only for x ^ X Q , i.e., outside the 
vortex cores. The field equation for ip derived from (Q) can be easily solved to yield 



(88) 



Upon taking the gradient of this equation, we obtain as velocity field v = (h/m)Vip 
the required form (46). Apart from a prefactor this expression coincides with the one 
for a cxt in (|87|). This is because a cxt can be written as a gradient of a scalar function 
[see (p6[)], When substituting the solution ( |88| ) back into the Lagrangian ( |84|) , which is 
tantamount to integrating out the field <p, we recover precisely the action (j48[). We thus 
see that the theory (^) with its linear coupling to a Chern-Simons field correctly re- 
produces the two-dimensional action of point vortices in an incompressible superfluid. 

Let us continue by considering the lowest order elastic scattering amplitude of two 
phonons calculated from the effective theory. In the frame where the sum of the mo- 
menta of the two incoming phonons is zero, the class of diagrams involving the ex- 
change of Chern-Simons quanta represents the scattering of a phonon from a vortex. 
An analogous situation arises in the context of Aharonov-Bohm scattering, i.e., elas- 
tic scattering of a charged nonrelativistic particle from an infinitely thin magnetic flux 
tube. It was pointed out by Bergman and Lozano [39 1 that such a scattering process can 
also be described by a quantum field theory consisting of a nonrelativistic |i/;| 4 -theory 
(with zero chemical potential) coupled to a Chern-Simons term — albeit minimally. 

At small energies and momenta, we can neglect the higher-order terms in the La- 
grangian (77) and restrict ourselves to terms at most quadratic in the field tp: 



£(2) 




l -(nd 0V ) 2 - c ^{hv v f + c 

-—doipVip • a I + £ C s- 
mc 



hd ip+ — {hV<p) 



2m 



(89) 



If we again introduce external vortices by replacing the Chern-Simons field a with (p7[), 
the field equation for ip(t, x ^ X Q ) becomes 



E%<p - c 2 \7 2 ip = 2—duVip ■ a c 



mc 



(90) 
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where, in accordance with the derivation of the hydrodynamic equation ([74]), we omit- 
ted a term proportional to do& cxt . Given the form ( p7| ) of a cxt , we recognize in the 
field equation ( p0[ ) derived from our theory, the hydrodynamic equation ([74]) with the 
contribution from the vortex motion ignored. Since Eq. ([74]) was used as starting point 
for the quantum-mechanical description of the elastic scattering of a sound wave from 
a vortex, we expect the effective theory (^9|) to be the appropriate starting point for the 
quantum-field-theoretic description. 

The propagators and vertices needed to calculate the scattering amplitude are read- 
ily obtained from that Lagrangian (|9|). In the gauge V • a = 0, the nonzero components 
of the Chern-Simons propagator follow as 

k k- 
waaAaaaaaa, : it& i0 (k) = -iW i(k)=-he ij j±, (91) 

where k is the wave vector. Note that there is no frequency dependence here. The 
phonon propagator reads 



P 



ih 



> : ihG(p) = —- 2 2 . , (92) 

n Pq — crp z + it] 

where r\ is a small positive constant that has to be taken to zero after the integration 
over the loop energy po has been carried out. In (|2|), p denotes the momentum. As 
vertices we read off 



% % GTl 

-rTifaq) = --^^(Poqi+PiQo) (93) 




% % GTl 

-r (p )(7 ) = --^ r ^p-q (94) 
h h m c 



P o 



1 en 
h mc 2 



-www : -— Po- (95) 



In [h3[|, we calculated the scattering amplitude A(9, |p|) to lowest order, with the result 



P 1 
-»* * *- 



^JP|) = -|44?^ (96) 
2 m i c L 1 — cos 
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where, as in (|7q), 9 is the scattering angle, and |p| the momentum of the incoming and 
scattered phonons. In deriving this, use is made of the following relations between the 
energy-momenta of the incoming and outgoing phonons 

Po=Po = qo = qo, |p| = |p| = |q| = |q| =Po/c, p = -p, q=-q, (97) 
and 

p-q = p 2 cos(#), (98) 

in the frame defined by the condition that the sum of the momenta of the incoming 
phonons be zero. The normalization £ of the external lines was determined in [|l3|] to 
be given by the dimensionless factor 



2ponA 



(99) 



where A is the surface area of the system. Apart from a kinematic factor, which is 
due to a different definition of the scattering amplitude in quantum mechanics and in 
quantum field theory, Eq. ( ^6|) agrees with the result (76) of Pitaevskii Jl6|]. Hence, 



using the effective theory (|89|), we have reproduced the Born approximation of the 
scattering of a sound wave from a vortex. 

For identical particles, one must add to the diagram in (^) the crossed diagram 
with the two outgoing lines exchanged. This yields the same result as in ( p6| ) with 
9 — ► 9 — ir. Adding the two contributions, we find 

A^(e,\p\)=i— cot(0), (100) 

777 

where we introduced the dimensionless parameter 

a = e|k|, (101) 



which accordi ng t o (|35|) was assumed to be small. Remember that p = Kk. The 
ratio h/rn in ( |100| ) is the circulation quantum. As observed by Sonin the result 
( |100| ) is identical to the lowest-order Aharonov-Bohm scattering amplitude, where the 
parameter a denotes the magnetic flux through the flux tube measured in units of $o> 
and m the mass of the scattered particle. The ratio h/m should now be interpreted not 
as the circulation quantum, but as the dispersion constant [ |40| ] which determines the 
dispersion of a matter wave with wave vector k describing a nonrelativistic particle of 
mass 777, 

w(k) = -^k 2 . (102) 

Z777 

In the quantum-field-theoretic approach to Aharonov-Bohm scattering [B9], the contri- 



bution (10C), which diverges for both 9 = and tt, is due to the exchange of a single 
Chern-Simons quantum. Since the Aharonov-Bohm scattering amplitude is known ex- 
actly, the next term in a perturbative calculation of this amplitude is known to be of 
order a 3 . 
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1 2 l-l I 



Figure 3: Contribution to the elastic scattering amplitude owing to the exchange of 
I Chern-Simons quanta. The incoming and outgoing phonons have an energy and 
momentum of the order cP and P, respectively. 



5 One-Loop Corrections 



In this section, we consider one-loop corrections to the tree amplitude (100). We are 



interested in the next order in an expansion in powers of the dimensionless quantity 



a introduced in (101). Because of the analogy with Aharonov-Bohm scattering, we 
expect this contribution to be of order a 3 and not of order a 2 . We calculate these 
contributions again using the effective Lagrangian (|8^). 

To identify the relevant Feynman graphs, we apply — as is commonly done in the 
context of effective field theories — dimensional analysis to the higher-order contribu- 
tions to the scattering amplitude. Following Weinberg we assume that the energies 
and momenta of the two incoming and the two outgoing phonons are small and all of 
the same order, po = cP, |p| = P say. As argued in that reference, the loop integrals 
will be dominated by contributions of the order P, so that one can make a perturbation 
in powers of P. More specifically, each energy and momentum appearing in the in- 
teraction vertices contributes a factor P. An internal phonon line therefore contributes 
a factor l/P 2 [see Eq. (p2|)], while an internal Chern-Simons line contributes a fac- 
tor 1/P [see Eq. |Tl])]. An integration over a loop frequency and loop wave vectors 
contributes a factor P 3 in two space dimensions. 

It is now easily shown that the exchange of I Chern-Simons quanta yields a con- 
tribution of the order of p 21 ^ 1 to the scattering amplitude (see Fig. |]), where we also 
took into account the energy-dependence contained in the normalization factor (^) of 
the external phonon lines. For the exchange of a single Chern-Simons quantum (/ = 1), 
we recover the momentum dependence given in ( |100| ). Since the two vertices involved 
are both proportional to the charge e with which the phonons couple to the Chern- 
Simons field describing the vortices, these diagrams are proportional to e 21 . Including 
also the other constants contained in the vertices and the propagators, we find that the 
dimensionless expansion parameter for this type of diagrams is a 2 . These diagrams 
represent the scattering of a phonon from a vortex, and are the contributions we wish 
to evaluate. 

For comparison, let us consider the case where instead of Chern-Simons quanta, 
phonons are exchanged. The interaction vertices involved in these contributions arise 
from the nonlinear terms in the effective theory which we have ignored in (|9|). Using 
similar arguments as the ones above, one can readily show that the exchange of I' 
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phonons yields a contribution to the phonon-phonon scattering amplitude of the order 



of P 



31' -1 



To be more specific, the dimensionless expansion parameter is here 



n 



(103) 



It is gratifying to see that the contributions which do not involve the Chern-Simons 
field describing the vortices have a different P-dependence. 

From the above dimensional analysis, it follows that to order a 3 , the contributions 
to the scattering amplitude are of the form 



A$ t (e,\p\)=i— [Co^+d^lpDa 2 ], 



ha 



(104) 



where naively the coefficient of the second term, C\ , is expected to be a function of the 
scattering angle 9 only, in the same way as the coefficient of the first term, 



C o (0) =cot(0), 



(105) 



depends only on 6. However, a one-loop diagram produces a nonanalytic term of the 
form (ha 3 /m) ln(|p|), so that C\ depends also on |p|. Whether such logarithmic terms 
will appear can a priori not be determined solely by dimensional arguments because 
that function is dimensionless. 

The calculation of the relevant one-loop corrections to the elastic phonon-phonon 
scattering amplitude is somewhat technical and therefore relegated to the Appendix. 
The final result for the coefficient C\{6, |p|) introduced in (104), which we split as 



Ci(0,|p|) :=Ci(d) + C„a(|p| 



(106) 



reads 



Ci{6) = -— 5 1 -+sm(i0)cos(i0)[±7r-2sm(|0)cos(±( 



+ sin J (±(9) 
+ cos 3 (i( 



1 - sin(±0 
cos(i0) 
sm(|0) 

1 + cos(i( 



1 / 6-TT 

~ 2 ^ + cos(i0) 

1 ( o 

2 V an(ie) 



(107) 



and 



I6ir 3 



dy- 



(108) 



This last contribution, which is obtained from Eq. ( A. 25 ) with the crossed term in- 
cluded, is logarithmically diverging. It is independent of the scattering angle. A closer 
inspection reveals that the divergence is an infrared divergence arising when the mo- 
mentum of one of the Chern-Simons propagators in the diagram depicted in Fig. ^ 
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Figure 4: Graphical representation of the O (a)-contribution Cq(9) and the O (a 3 ) 
contribution C\ (9) to the elastic phonon-phonon scattering amplitude. 



tends to zero. To regularize the integral, we give the Chern-Simons field a small mass 
A*; Cna(|p|) then becomes 



C„a(|p|) = - 



1 



dy- 



32^ J " („ + !)(„_ _ ^2/^)2 

i f°°. y-i 



32^^ *(y + l )(2/ _^l_ M 2 c 2 /p2) 2 



1 



16tt 3 



In 



(109) 



This introduces an arbitrary scale factor into the problem. Frequently, such a scale 
factor arises instead when ultraviolet divergences are present. Logarithmic divergences 
of this kind have to be eliminated by a redefinition of certain parameters of the effective 
theory at the expense of the appearance of an arbitrary renormalization scale . The 
nonrelativistic effective theory considered here is ultraviolet finite to this order, and 
needs, therefore, not to be renormalized. The arbitrary scale factor /i arises here from 
the infrared region. 

A graphical representation of the analytic part of the one-loop result, C\(9), to- 
gether with the tree result is given in Fig. ^. It shows that whereas the tree contribution 
Cq(8) is antisymmetric about 9 = ir/2, the O (a 3 ) -contributions are symmetric. Both 
are seen to diverge for 9 — and ir. 



6 Discussion 

We have extended the analysis of the quantum field theory recently proposed by these 
authors [[□]] to describe vortices in a superfluid film at the absolute zero of temperature. 
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Figure 5: One-loop diagram I contributing to the scattering amplitude. 

The theory consists of the effective action of phonons coupled to a Chern-Simons term. 
Two of its salient features are that, first, the effective action describing the phonons is 
invariant under Galilei transformations — as is required for a nonrelativistic system — 
and, second, the coupling to the Chern-Simons terms is linear and not minimal. We 
have demonstrated that various known facts about two-dimensional quantum vortices 
are correctly reproduced by the quantum field theory. 

We have further shown how it can be used to calculate one-loop contributions to 
the amplitude for the elastic scattering of phonons from a vortex. We have applied di- 
mensional analysis of the type commonly used in the context of effective field theories 
to identify the relevant Feynman graphs and have calculated these. We have shown 
that the one-loop corrections are ultraviolet finite, and that an arbitrary scale factor was 
introduced by an infrared divergence. 

By minimally coupling the Goldstone field ip to electrodynamics, the theory dis- 
cussed here may be extended to describe vortices in a superconducting film. 
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Appendix 

In this Appendix, we calculate the contributions to the elastic phonon-phonon scatter- 
ing amplitude proportional to a 3 . These contributions arise from the three one-loop 
diagrams depicted in Figs. |5]-[7|. Since the dimensional analysis of the previous section 
already fixed the dimensionful prefactors appearing in the scattering amplitude, we can 
simplify matters by setting all dimensionful parameters to one, 

h = c = m = fi = l. (A.l) 

The Feynman rules (^T|)-(p4|) applied to the diagram depicted in Fig. ||, which we 
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refer to as I, yield 



i = C 4 



^G(r)a(r)a il AL-)a r ^k)T ll{l ,-r)T r (iK - r)T,jr. -qKM: ■<}). 



(A.2) 

where because of conservation of energy and momentum at the vertices we have in 
addition to <^T^ the following relations 

r a = 2p ~ r , k = p - r , k = r -p , f = -r, k = p r, k = r q. (A.3) 



The summation over the indices in (A.2) (using an Euclidean metric) gives 



G t _ L y{k)G p(T (k)T^{p, -r)T u {p, -r)T p (r, -g)r CT (f, -q) = 

[p • r(p fxk-f pxk)+p- f (po r x k + r p x k)] 



k 2 k 2 



x [q • r(p fxk-f qxk) + q- f (p r x k + r q x k)]. (A.4) 



Substituting this in (A.2) and using the relations (A.3), we obtain the expression 



1 = 4 



d-V 



prpxrqrqxr 



(2?r) 3 (rl - r 2 + irf)[{2p Q - r ) 2 - r 2 + irj\ (p - r) 2 (q - r) 2 



(A.5) 



With the help of contour integration, the integral over the loop energy tq is readily 
carried out with the result 



d 2 r 



prpxrqrqxr 



(2^) 2 |r|(r 2 -p 2 )(p-r) 2 (q-r) 



(A.6) 



The integral over the angle a is most easily carried out by introducing the variable 
z = exp(ia) and performing a contour integration along the unit circle. The poles one 
encounters are located at 



z = 0, z = y, z = 1/y, z = e l6 y, z = e iS /y, 



where y stands for 



y 



Po 



(A.7) 



(A.8) 



Since only those poles lying inside the unit circle are to be included, the regions < 
y < 1 and 1 < y < oo have to be treated separately. Adding the two contributions, we 
find 



I 



s(26>) + 2cos(26V - [1 + 2cos(6)]y 4 
1 - 2cos(% 2 +y 4 



The remaining integral over y is elementary and yields 

-ipl ( 10 14 



I 



16tt 



■cos(0) +4cos 3 (i(9) ln[tan(i6»)] 



(A.9) 



(A. 10) 
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Figure 6: One-loop diagram II contributing to the scattering amplitude. 



This result is valid for < 9 < 2tt. 

We continue with the evaluation of the second one-loop diagram depicted in Fig. ||. 
For this diagram II we find 



II = c 4 



G(r)G(?)G^(k)G pa Ck)T P ( P , -r)r„(f, -q)T p (p, -r)r a {r, -q), 

(A.ll) 



where now 

?~o = 7"o, k = p -r , k = p -r , f = r-p-q, k = p r, k = q r. (A.12) 
Explicitly, 

G llv (k)G pa (k)T fl (p ) -r)T„(r, ~q)T p (p, ~f)T a (r, -q) = 

[p • r(p Q fxk-f qxk) + q- f (p r x k + r p x k)] 



k 2 k 2 



x [p • f (p rxk + r qxk) + q- r(p f xk-f px k)]. (A.13) 



Substituting this in ( |A.l 1[ ) and using the relations (A.12), we arrive at 

(Po + ro) 2 1 



II 



1 f d 3 r 

Ap 2 Q J (2?r) 3 (rl - r 2 + ^)[r 2 - (r - p - q) 2 + irj\ (p - r) 2 (q - r) 2 
x [p • r(p xq + qxr) q ■ (r — p — q)p x r] 
x [— q • r(— p xq + pxr)+p-(r — p — q)q x r]. (A.14) 



Let us first carry out the integral over the loop energy tq: 



dr 



(Po + r ) 2 



i (po 



i (po - W\ 



2tt (tq — r 2 + «7?)(tq — v 2 + if]) 



2 r r 2 - v 2 



2 v v 2 



r 2 ) 



(A.15) 



where we introduced the abbreviation v = r p q. The first term in (A.15) gives 



a contribution Hi to II and the second term a contribution II2, such that II = Hi + II2. 
It turns out that both contributions are identical, so that II = 2IIi . We perform the 
angle integration again by introducing the variable z = cxp(ia) and carrying out a 
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Figure 7: One-loop diagram III contributing to the scattering amplitude. 
contour integration along the unit circle. Besides poles located at z = y, e l8 y and 



z 



1/y, e /y, with y defined in (A. 8), there are also two poles at 



l+e l6 ± ^( e ^ + l) 2 - 4e l V 
21,2 = ^ • (A.16) 

The contributions of the first two pairs cancel, so that we only have to consider the 
contributions stemming from the poles z\ and For y > cos(^8) these poles lie 
precisely on the unit circle, \zi^\ = 1 and it is a priori not clear how to proceed. A 
closer (numerical) inspection reveals that the contributions from z\ and z 2 cancel there. 
This leaves us with region y < cos(i#), where \zi\ < 1 and \z 2 \ > 1. In this way, we 
obtain with t = cos(^8) the result 

n _ -ipl f l Ay U{l-y)[yHe-l)+e] 



16tt 



16tt Jo (l + y)y/t 2 -y 2 

{-4sin 2 (±6>) cos 2 (i6») + | cos( ±6>)[4 - 5 cos 2 (±6>) - cos 4 (±(9)] 
+20cot(i0)cos(0)}, (A.17) 

which is valid for —ir<9<ir. 

The last one-loop contribution to the phonon-phonon scattering amplitude is de- 
picted in Fig. 0. It stands for 

HI = C 4 J ^G{r)G{f)G^{k)G pa {k)Y^p, -r)T v (-q, -r)T p (p, r)T a (~q, f), 

(A.18) 

where now we have the relations 

7~ = -r , k = p -r , k = p + r , r = p-q-r, k = p-r, k = r-p. (A.19) 
More explicitly, 

G pv (k)G pa (tyT^p, -r)Y v (-q, -f)T p (p, r)T a (-q, f) = 

— ~- [p • r(p fxk + r qxk)-q - f (p r x k + r p x k)] 
k 2 k 2 

x [p • r(p fxk — f qxk) — q - f (p r x k - r p x k)]. (A.20) 



26 



In this way, ( A..18) becomes 
III 



d 3 r 



-1 



Ap 2 J (27r)3(rg 



2 9 

Po -»o 



*»7)[»o _ (P- 



irj\ 



(A.21) 



[p • r(p xq + qxr) + q(p q r)p x r] 5 



(p-r)4 

We again first carry out the integral over the loop energy j-q, 



dr Pq ~ r o * 1-2 ^Po 



vl 



2tt (rg 



i?7)(r 2 



2 Irlfr 2 



2 Iwlfw 2 



(A.22) 



where w = r + q p. The first (second) term gives a contribution IIIi(2), such that 
III = IIIi + III2. The angle integral is a gain evaluated by a contour integration along 
the unit circle. With y as defined in (A. 8), the integrand of IIIi has poles at z = y, 1/y 
as well as at 



Z3A 



1 - e %e T V(e< 9 - 1) 



0,,,2 



2y 



while the integrand of III2 has poles at z = e l6 y, e 10 /y and at 

-1 



1 



Ae l0 y 



e, .2 



2>J 



(A.23) 



(A.24) 



In III2, we made the change of integration variable r — > w. The contributions from 
the poles z — y,l/y, e l9 y, e 10 jy can be calculated in a straightforward manner. They 
cause a logarithmic-diverging integral 



IHn 



^sin 2 (^ 

47T V2 



dy- 



1 



y- 



1 



(A.25) 



which is discussed in the main text. The remaining poles z 3 , z 4 , z 5 , and z e which we 
have to consider lie precisely on the unit circle for y > sin(|#). A careful numerical 
analysis reveals that their contributions cancel here. This leaves us with the region 
y < sm.(^6), where only the poles 2:3 and zq have to be included since |z3,e| < 1 while 
1^4 5] > 1. Both Z3 and zq give the same real part, while their imaginary parts cancel. 
With s denoting sin(i#), we finally obtain for diagram III 



in - m na 
i 

16tt 



(A.26) 



4s( 2 / 2 -l)[-s 2 + 2 / 2 + (l-s 2 )y 4 



V s " 2 - v 2 



2sy 2 — \J s 1 — y 2 ^j + (^s — yj s 2 — y 2 ^j 



x 8s 4 



isV + y* 



8s 3 \/s 2 - y 2 + <isy 2 v / s 2 - y j = 



2 - 



{2n\ sin(±0)l cos(|0) - || sin(±0)l [4 - sin 2 (l0) - sin 4 (l0)] } 



where 8 is to be restricted to the values —tt < 8 < tt. 
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Since we are dealing with identical particles, we have to include the crossed di- 
agrams whose contributions are obtained from the uncrossed diagrams by replacing 
9 — > 9 — 7r mod(27r). More specifically, restricting ourselves to the values < 9 < it, 



we have to replace 9 with 9 + n in the result (A. 9) of diagram I, and with 9 — tt in the 
results (A. 17) and ( |A.26 ) of diagram II and III. The final expression for the one-loop 



contributions to the elastic phonon-phonon scattering proportional to a is given in 



(107) and (109) 
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